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Abstract 

In two-photon collisions at LEP2 and a future e+e^ linear collider heavy quarks (mainly 
charm) will be pair-produced rather copiously. The production via direct and resolved 
photons can be distinguished experimentally via a remnant-jet tag. We study correla- 
tions of the heavy quarks at next-to-leading order in QCD in the direct channel, which 
is free from phenomenological parton densities in the photon. These correlations are 
therefore directly calculable in perturbative QCD and provide a stringent test of the 
production mechanism. 
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1. Introduction 



The production of heavy quarks in two-photon colhsions has interesting aspects. Each 
of the photons can behave as either a pointhke or a hadronic particle Consequently 
one distinguishes in such collisions direct- (both photons are pointlike), single resolved- 
(one photon is pointlike, the other hadronlike), and double resolved (both are hadronlike) 
production channels. The resolved channels require the use of parton densities in the photon, 
whereas the production via the direct channel is free of such phenomenological inputs and 
depends only on the QCD coupling and the heavy quark mass. The heavy mass provides 
the hard scale for the perturbative analysis and ensures that the separation into direct 
and resolved production channels is unambiguous even at the next-to-leading order (NLO) 
level. Hence production via the direct channel is directly calculable in perturbative QCD 
(pQCD) and in principle the best way for examining the validity of such an analysis and 
for confronting the pQCD prediction with experiment. 

Two-photon collisions can be investigated at e^e~ colliders, where a large number of 
equivalent photons is generated. Charm quark production in two-photon collisions has been 
analysed in many experiments. One has mainly studied the reaction e'^e~ e^e~ D*^X 
with neither outgoing lepton tagged ( "no-tag" ) , because it proceeds predominantly via the 
fusion of two equivalent photons to produce open charm (77 — > cc). The existence of the 
D*^ has been inferred either from direct reconstruction Q or from unfolding the distribution 
of soft pions 1^] resulting from its decay. There have in addition been studies that use soft 
leptons [Q] and kaons |^] to tag charm quarks. 

Due to the low experimental acceptance of heavy quark production in two-photon col- 
lisions this reaction has been studied also theoretically at next-to-leading order in QCD 
only in the single-particle inclusive case. Ref.Q concentrated on the no-tag case, and ref.Q 
on the case where one of the outgoing leptons is tagged. At LEP2 and a future e^e~ lin- 
ear collider (NLC) the higher cms energy and large luminosity will lead to fairly copious 
production of charm quark pairs. Thus it will become possible to measure both heavy 
quarks and analyse their correlations. The study of these correlations constitutes a more 
comprehensive test of the theory and is our purpose in this letter. Heavy-quark correla- 
tions have been investigated theoretically also in hadroproduction Q, photoproduction |^ 
and electroproduction [^], and experimentally in We concentrate here on the no-tag 
case and, to eliminate the uncertainties related to the parton densities in the photon, on 
the direct channel only. Note that the TOPAZ collaboration Q has recently shown that 
the direct channel may be isolated experimentally from the resolved ones by detecting the 
photon-remnant jet, present in the resolved channels only. 

The paper is organized as follows: In section 2 we describe our method of calculation 
and in section 3 we show heavy quark correlations for LEP2, and a NLC at a center of mass 
energy of 500 GeV. We conclude in section 4. 
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2. Method 



In this section we describe the method we used to calculate the QCD corrections to the 
process 

lih) + 7{k2) ^ Q{pi) + Q{P2) , (1) 

where Q{Q) is a heavy (anti)-quark. We want to have full exclusive information about the 
final state. Our method is a special case of a more general method for performing exclusive 
higher order QCD calculations 

The Born process (|l|) is described by the differential cross section 



da(0) 



S \Ul ti tiUi 




P1-P2). (2) 



Here eg is the charge of the heavy quark in units of e, A'^c = 3 the number of colors, and m 
the mass of the heavy quark. The kinematic invariants are defined by 

s = {ki + k2f, ti = {ki-pif -m^, ui = {ki - p2f - . (3) 

The virtual QCD corrections to the Born process consist of the interference between the 
Born amplitude (depicted e.g. in Fig. Al in [^) and its one- loop corrections. Explicit results 



have already been presented in [13| and we will not repeat the details of the calculation here. 
We merely note that we regularized the ultraviolet (UV) and infrared (IR) singularities that 
occur in the virtual corrections by working in d = 4 — 2e dimensions, and absorbed the UV 
singularities via mass renormalization in the on-shell scheme. We are then left with only 
IR singularities, which appear as 1/e poles and factorize into a universal factor multiplying 
the Born differential cross section, eq.(||). 

The bremsstrahlung corrections at NLO are due to the radiation of a gluon from one of 
the heavy quarks 

l{ki) + 7(^2) ^ Q{pi) + Q{P2) + 9{ks) . (4) 

Since our method here is a little different from what was done previously in the literature, 
we give a few more details. Note first of all that when a gluon is radiated from the heavy 
quark, no collinear singularity occurs, because it is shielded by the heavy quark mass. We 
divide up the phase space into a "soft" region and a "hard" region. The soft region is 
defined by the condition 

< Si3, S23 < Smin (5) 

where Sjs = 2pi ■ k^ (i = 1,2) and Smin is an arbitrary cut-off, to be chosen small. The hard 
region is the complementary one. 

In the hard phase space region, one can work in 4 dimensions and perform the phase 
space integrations numerically, allowing for easy implementation of experimental cuts. As 
is well known, in the soft region both the phase space and the matrix element factorize in 
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the limit of small Smin. In both cases, one of the factors contains the quantum numbers 
of the gluon, and the other is only related to the lower order process. As a consequence, 
one may perform the integral of the momentum of the gluon in this region analytically in 
d-dimensions. Specifically one must do the integral 

^TTasCp j dPS(soft)i^(soft) . (6) 
with the color factor Cp = {N'^ — 1)/{2N). Here the soft gluon phase space factor is 

"^^^^'"^^^ = leJra- g) d^i3dg23(^/3)^^-H^i2^i3^23 " m^sj, + 4)]-^ (7) 

where /3 = ^/T^^AmFJs and S12 = 2pi ■ p2. Note that the expression in square brackets 
must be positive. The soft gluon matrix element factor can easily be found in the eikonal 
approximation, and is 



i^(soft) = 4 ( ^1^^13^^3-r.^(^f3 + ^i3) \ _ 
\ '5l3'^23 / 



Thus, upon combining both factors and integrating over the range (^), one obtains a univer- 
sal factor multiplying the differential Born cross section (^). This factor contains 1/e poles 
which cancel against those originating in the virtual corrections. The soft contribution to 
the fully differential cross section can then finally be written as 



where 



dCT(^)(soft) = SF(s,m2,Smin) dcj(°) (9) 



-2(ln(l - x) +ln(l +x) - Inx) + 1-/3 

lnx(l + 21nil±^Kl^) (10) 

Li2 ^) -Li2(l-x2; 

m^, f(l + x)(l-x) 3s ( 2m2\, ^ 
+— In X ^ '-^ + 1 In X 




Here x = (1 — /?)/(! + /3) and Li2(2;) is the dilogarithmic function as defined in |14|. 

Finally, one is left with a two-to-two particle contribution (consisting of the Born and 
soft-plus-virtual corrections) and the two-to-three particle contribution in the hard region. 
Each contribution depends on the theoretical cut-off Smim but as long as Smin is small enough 
compared to the typical scale of our process, the sum does not. This we checked explicitly. 
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3. Results 



Using the method described in the previous section we have constructed a Monte Carlo 
program for the reaction 77 QQ for direct photons, including the complete 0{as) correc- 
tions, which is fully exclusive in all final state particles. We checked that we could reproduce 
the results in for the total cross section and single particle transverse momentum (pt) and 
rapidity (y) distributions for the direct channel. We only present results for charm quark 
production because the bottom quark production rate is very much reduced in two-photon 
collisions due to charge and phase space suppression. 

We first list the default choices we made for various parameters for producing the results 
shown in the rest of this section. To compute we used the two loop expression with 
Aqqj-, = 0.215 GeV and nu active flavors, where nu is the number of flavors with mass less 
than the renormalization scale. For the charm quark mass we used 1.5 GeV. The center 
of mass energy was chosen to be 175 (500) GeV for LEP2 (NLC). For the renormalization 



scale we took /i = ym^ + (PtiQ) +Pt(Q))/2. In the present process the choice of scale 
only affects the value of a^. We used the Weizsacker- Williams density of ||l^ with an anti- 
tag angle ^max of 30 (175) mrad for the case of LEP2 (NLC). At the NLC beamstrahlung 
is expected to play an important role, so we include its effect here by adopting for its 
spectrum the expression given in [16|, with parameters Tgfr = 0.039 and cr^ = 0.5 mm [17| 
corresponding to the TESLA design. For the NLC we will as default coherently superimpose 
the Weizsacker- Williams density and the beamstrahlung density, in order to incorporate the 
case where one photon is of beam- and the other of bremsstrahlung origin. 

For most results we have not used charm-to-D meson fragmentation function. For the 
cases that we do, which we indicate explicitly, we employed the Peterson et al. parametriza- 
tion |l|] 

D{z) = ^ ^ (11) 

zil — 1/z — e/(l — z)Y 

with e = 0.06 the value given in [^] for the case of charm. Our interest when including 
the fragmentation function lies mainly in how it changes the shapes of distributions, rather 
than their normalization. Hence we choose such that /q^ dz D{z) = 1. 

We will only present one single particle distribution here, since such distributions have 
already been studied in [I52C]. Fig.l shows the single particle pt distribution at LO, NLO, 
and at NLO with fragmentation. We see that inclusion of NLO corrections decreases the 
cross section at large pt and enhances it at small pt, and that the application of the frag- 
mentation function softens it considerably. 

Turning to correlations, we begin by showing distributions which allow a comparison 
between the LO and NLO calculations. In Fig. 2 we show the cross section versus the 
invariant mass Mqq of the heavy quark pair for LEP2 and NLC at both LO and NLO. 
Notice in Fig. 2 the sizable difference that occurs at both small and large invariant masses 
when including the NLO corrections. This can be understood as follows. Consider flrst 
the situation where the two photons collide with all the momentum of their parent leptons. 
Denoting the invariant mass of the heavy quark pair in this case by Mqq , then at LO 
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Figure 1: Single charm quark pt spectrum at LEP2, comparison of LO, NLO and NLO 
with fragmentation. 




Figure 2: Mqq distribution for charm for LEP2 and NLC at both LO and NLO. 



Mqq is fixed at ^/s^. At NLO it may assume smaller values, and there the cross section 
is positive. For Mqq = y/sT^ one has at NLO also a negative contribution coming from 
the virtual graphs. To go back to the case of LEP2 and NLC we must fold with the 
photon spectrum. A given Mqq value then contributes to the spectrum for Mqq under the 
restriction Mqq < Mqq, so that at large Mqq the LO spectrum is mainly modified by the 
negative contribution at Mqq = s/Sjj, and at small Mqq by the positive contributions at 
smaller Mqq. 

In Fig. 3 wc show the AR distribution, defined by AR = \/ (A(/>)'-^ + {Ar])'^, at LO and 
NLO for both LEP2 and the NLC. Here is the azimuthal angle between the charm and 




Figure 3: AR distribution for charm and anti-charm quark at LEP2 and NLC. 

anticharm in the plane transverse to the beam axis and Ar] is the pseudo-rapidity difference 
of the two heavy quarks. At LO AR > it, but at NLO AR may also assume values below 
that. Note that NLO effects seem to be mostly active for AR <, 4. 

We now show two distributions which are only non-trivial at NLO (and higher orders). 
In Fig. 4 we present the pt distribution of the charm-anticharm pair, and in Fig. 5 the az- 
imuthal correlation between the two heavy quarks. We also show in Fig. 4 the NLC curves 
with only beamstrahlung photons and with only Weizsacker- Williams photons for the pur- 
pose of comparison. One observes in Fig.4 that at the NLC charm pairs produced by 
beamstrahlung photons prefer to have a lower pt than those due to WW equivalent pho- 
tons. This is a consequence of the TESLA beamstrahlung spectrum, which is enhanced 
at small z and depleted at large z compared to the WW spectrum (2 is the momentum 
fraction of the photon relative to its parent lepton). 
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Figure 4: pt{cc) distribution for charm and anti- charm quark at LEP2 and NLC. 
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Figure 5: Acf) distribution for charm and anti-charm quark at LEP2 and NLC. 



In Fig.5 we see that the distributions are all quite uniform. We observe however 
that at the NLC for the case of charm the beamstrahlung contribution dominates the WW 
one. 

Finally we comment on the consequences of choosing different values of the renormal- 
ization scale /x and the charm mass m. To see how Figs. 1-3 change when varying /x one can 

simply rescale the differences between the LO and NLO curves according the change in as, 
whereas in Figs. 4 and 5 the whole curve will change by an overall factor. We further remark 
that choosing a different value for m changes mainly the normalizations of the curves shown 
in this section, but not their shapes. 

4. Conclusions 

In this paper we have presented a NLO calculation of heavy quark production in direct 
two-photon collisions. We have described our calculation method and presented numerical 
studies of various correlations between the heavy quarks. We observed that the inclusion of 
the NLO corrections significantly modifies the shapes and normalizations of the distributions 
we studied. Experimentally such studies will be challenging at LEP2 due to the low charm 
acceptance, but they are certainly feasible at a future e'^e" linear collider. 

References 

[1] T.F. Walsh and P.M. Zerwas, Phys. Lett. B44 (1973) 195; E. Witten, Nucl. Phys. 
B120 (1977) 189. 

[2] W. Bartel et al. (JADE), Phys. Lett. B184 (1987) 288; M. Alston-Garnjost et al. 

(TPC/27), Phys. Lett. B252 (1990) 499; W. Braunschweig et al. (TASSO), Z. Phys. 
C47 (1990) 499; R. Enomoto et al. (TOPAZ), Phys. Rev. D50 (1994) 1879; D. Buskulic 
ct al. (ALEPH), Phys. Lett. B255 (1995) 595; T. Nozaki (AMY collab.), Proceedings 
of Photon'95, Sheffield (1995), S. Cartwright and D. Miller cds. (World Scientific). 

[3] R. Enomoto et al. (TOPAZ), Phys. Lett. B328 (1994) 535. 

[4] S. Uehara et al. (VENUS), Z. Phys. C63 (1994) 213. 

[5] R. Enomoto et al. (TOPAZ), Phys. Lett. B341 (1994) 99. 

[6] M. Drees, M. Kramer, J. Zunft and P.M. Zerwas, Phys. Lett. B306 (1993) 371. 

[7] E. Laenen, S. Riemersma, J. Smith and W.L. van Neerven, Phys. Rev. D49 (1994) 
5753; E. Laenen and S. Riemersma, in preparation. 

[8] M. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B373 (1992) 295. 

[9] S. Frixione, M. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B412 (1994) 225. 

[10] B. Harris and J. Smith, Phys. Lett. B353 (1995) 535; Nucl. Phys. B452 (1995) 109. 



8 



[11] M. Aguilar-Benitez et al. (LEBC-EHS coll.), Phys. Lett. B164 (1985) 404 and Z. 
Phys. C40 (1988) 321; M.I. Adamovich et al. (Photon Emulsion coll.), Phys. Lett. 
B187 (1987) 437; S. Aoki et al. (WA75 coll.), Phys. Lett. B209 (1988) 113; S. Barlag 
et al. (ACCMOR coll.), Phys. Lett. B257 (1991) 519; K. Kondama ct al. (E653 coll.), 
Phys. Lett. B263 (1991) 579; M.P. Alvarez et al. (NA14/2 coll.), Phys. Lett. B278 
(1992) 385; V. Arena et al. (E687 coh.), Phys. Lett. B308 (1993) 194. 

[12] W.T. Giele and E.W.N. Glover, Phys. Rev. D46 (1992) 1980; W.T. Giele, E.W.N. 
Glover and D.A. Kosower, Nucl. Phys. B403 (1993) 633; S. Keller and E. Laenen, in 
preparation. 

[13] W. Beenakker, H. Kuijf, J. Smith and W.L. van Neerven, Phys. Rev. D40 (1989) 54. 
[14] L. Lewin, Polylogarithms and Associated Functions (North-Holland) Amsterdam, 1983. 
[15] S. Frixione, M. Mangano, P. Nason and G. Ridolfi, Phys. Lett. B319 (1993) 339. 
[16] P. Chen, T. Barklow and M.E. Peskin, Phys. Rev. D49 (1994) 3209. 
[17] D. Schulte, private communication. 

[18] C. Peterson, D. Schlatter, I. Schmitt and P. Zerwas, Phys. Rev. D27 (1983) 105. 
[19] J. Chrin, Z. Phys. C36 (1987) 163. 

[20] M. Cacciari, M. Greco, B.A. Kniehl, M. Kramer, G. Kramer and M. Spira, preprint 
DESY 95-205. 



9 



